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ILLUSTRATIONS 


Figure 

1  Geometry  of  stellatron  accelerator  showing  conducting 
toroidal  cavity  and  externally  applied  magnetic  fields. 

2  Frequencies  and  growth  rates,  denoted  by  wr  and  r, 
respectively,  obtained  from  the  dispersion  relation, 

Eq.  (9),  for  Table  I  parameters  with  b  =  8.8,  I  =  10  kA 

m  =  22. 

3  Growth  rates  of  instability  versus  quadrupole  mode- 
number  m  (which  must  be  an  even  integer)  for  parameters 
in  Table  I,  and  b  =  8.8  cm,  I  =  10  kA. 

4  Energy  in  n  =  22  field  components  (n  =  toroidal  mode- 
number)  versus  time  showing  development  of  i  -  0,  m  =  2 
instability  in  Figure  2. 

5  Contour  plot  of  n  =  0  component  of  toroidal  self- 
magnetic  field  B  during  early  nonlinear  growth 

of  an  i  =  0,  m  =20  instability.  The  contours  have 
linearly  increasing  values  from  A  to  G. 
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Growth  rate  of  instability  versus  beam  current  for 
Table  I  parameters  with  m  =  22,  b  =  9.0. 


I.  INTRODUCTION 


Betatrons  and  other  recirculating  accelerator  designs  have  been  studied 
in  recent  years  as  compact  accelerators  for  high-current  electron  beams.1  6 
Conventional7  and  modified1  betatrons  require  precise  matching  of  the  beam 
energy  to  the  vertical  magnetic  field  if  a  fixed  major  radius  is  to  be 
maintained.  To  overcome  this  restriction,  strong  focusing  in  the  form  of  a 
helical  magnetic  quadrupole  has  been  added  to  the  modified  betatron 
configuration.3*8  This  considerably  increases  the  tolerance  of  the 
device  to  mismatch.  In  this  paper,  we  show  that  helical  quadrupole 
focusing  can  lead  to  an  electromagnetic  instability  arising  from  a 
three-wave  interaction  between  the  static  helical  field,  a  transverse  mode 
on  the  beam,  and  a  transverse-electric  (TE>  waveguide  mode.  This  behavior 
is  reminiscent  of  that  which  occurs  in  planar  geometry  in  the  presence  of  a 
rippled  magnetic  field.9  Unlike  the  f ree-electron  laser  instability,  which 
was  recently  analyzed  in  the  presence  of  helical  quadrupole  focusing,10 
longitudinal  bunching  of  the  beam  plays  no  significant  role  in  the  insta¬ 
bility  which  we  describe  here.  Also,  the  instability  is  essentially 
independent  of  the  radius  of  curvature  of  the  device.  It  is  thus  not 
related  to  the  negative-mass  instability,  for  example.  The  main  effect  of 
finite  radius  of  curvature  is  to  discretize  the  toroidal  mode-numbers. 

This  can  be  important,  since  if  the  instability  width  is  narrow  enough,  it 
can  fall  between  two  allowed  mode-numbers  and  disappear  (see  Section  II). 

In  Section  II,  we  give  an  analytic  theory  of  the  instability.  In 
Section  III,  we  obtain  a  simplified  dispersion  relation  which  yields  an 
instability  criterion.  Section  IV  details  comparisons  between  the  analytic 
theory  and  three-dimensional  particle  simulations.  Section  V  gives  our 
conclusions. 

Throughout  the  paper,  we  use  convenient  units  where  the  electron  charge 
e  and  mass  m,  and  the  velocity  of  light  c,  are  scaled  out.  Lengths  are 
normalized  to  c/u>g  =  1  cm,  frequencies  to  ujq,  velocities  to  c,  densities  to 
WQm/4Tre^,  electric  and  magnetic  fields  to  mcug/e. 


II.  ANALYTIC  MODEL  OF  INSTABILITY 


The  geometry  and  coordinate  system  we  use  are  illustrated  in  Figure  1. 
An  electron  ring  of  major  radius  R  and  minor  radius  a  circulates  around  the 
center  of  a  conducting  torus  whose  major  and  minor  radii  are  R  and  b, 
respectively.  The  externally  applied  magnetic  fields  consist  of  a  vertical 
field  Bz  «  r~s,  where  s  is  the  external  field  index  (we  assume  s  =  1/2) 
a  solenoidal  field  B0  ,  and  a  helical  quadrupole  field  Bq.  This  is  the 
configuration  of  the  stellatron  accelerator.3  A  configuration  where  the 
helical  quadrupole  field  is  generated  by  two  current-carrying  wires  instead 
of  four  also  has  been  proposed.6  We  assume  that  the  beam  can  be  modeled  as 
a  string  of  rigid  disks.  This  means  that  we  follow  the  transverse  motion 
of  the  beam  centroid,  and  ignore  any  internal  dynamics.  This  treatment  is 
valid  provided  the  drift-tube  minor  radius  is  much  greater  than  the  beam 
radius,  b  »  a,  since  the  fields  we  need  to  consider  have  a  transverse  scale 
length  -b.  We  ignore  the  perturbed  e-motion  of  the  beam  under  the 
assumption  of  relati visitic  stiffness  in  this  direction,  i.e.,  y  >  1,  where 
y  is  the  beam  relativistic  factor. 

In  equilibrium,  the  matched  value  of  the  vertical  magnetic  field  is 

=  yV./R,  where  V.  is  the  toroidal  beam  velocity.  Linearizing  about  the 
z  w  o 

equilibrium  position,  we  obtain  the  following  equations  of  motion  for  the 
perturbed  coordinates  r,  z  of  the  beam  centroid.3 

•  1  O*  O  M  «•  *  * 

r  +  Ynzr  "  fiez  +  ufiz(rcos  17,6  +  zsin  me)  =  -Er  ■  veBz  ♦  (l) 

z  +  Yfizz  +  fier  +  unz(rsin  1,16  -  zcos  m0)  =  -EZ  +  veBr  »  (2) 

where  a  dot  denotes  the  total  time  derivative  3/9t  +  (V0/R)9/9e,  nz= 

B  /Y,  =  B  /y ,  u  denotes  the  amplitude  of  the  quadrupole  field  index,  and 

_  Z  m  u  0 

E^,  B  denote  the  perturbed  self  electric  and  magnetic  fields  acting  on 
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Figure  1.  Geometry  of  stellatron  accelerator  showing  conducting 
toroidal  cavity  and  externally  applied  magnetic  fields. 


the  beam  at  the  perturbed  position.  Letting  c  =  z  +  ir,  Eqs.  (1),  (2)  can 
be  written 


12  . 

^ + 


_2  -im6  r“ 

vQ  e  5*  =  -[E 


+  iE  + 
z  r 


ive(Bz  + 


1Br)] 


(3) 


where  *  denotes  the  complex  conjugate.  To  calculate  l,  J5,  we  use  a  Green's 
function  approach.  Maxwell's  equations  can  be  written 


v  x  v  x  E 


u»2E 


iwJ 


> 


(4) 


£  *  .  (5) 

where  we  assume  an  e  1wt+1te  dependence  for  the  perturbed  fields  and  cur¬ 
rents.  The  solution  to  Eq.  (4)  can  be  written  down  using  a  dyadic  Green's 
function,11  constructed  from  the  solutions  to  the  homogeneous  counterpart 
of  Eq.  (4).  The  expressions  obtained  are  infinite  series.12  The  expres¬ 
sion  for  Er,  for  example,  evaluated  at  the  location  of  the  beam,  is 


E  =  -va)(o)  -  jw  )r 

r  v  zv 


I 

n 


I 

n 

1 — T 

/nL 


n  +  i  /R 


~7 

OJ 


- 1 
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*2/R2  I 


(X2  +  *2/R2  -  w2)(x2  +  tZ/R2) 


8E°  - 


+  0 


where 


(6) 


J2(nb) 


b2 

T 


[J'(xb)]' 
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w  ere  is  a  Bessel  fui,  :tion,  and  where  n  denotes  the  roots  of  JJ(nb)  =  0, 
X  denotes  the  roots  of  J^(Xb)  =  0,  and  0(b/R)  denotes  toroidal  corrections, 
which  we  neglect  under  the  assumption  R  »  b.  Also,  v  is  Budker's  param¬ 
eter,  i.e.,  the  number  of  electrons  per  unit  length  of  the  beam  times  the 
classical  electron  radius.  The  last  term  in  Eq.  (6)  results  from  evalu¬ 
ating  the  equilibrium  radial  electric  field  E°  at  the  perturbed  position. 


In  deriving  Eq.  (6),  we  have  used  the  rigid-disk  approximation  to  write 
jp  =  _pr  =  _  tn^jpr,  where  p  is  the  beam  density.  Similar  expressions 

are  obtained  for  the  other  fields.  To  simplify  the  analysis,  we  extract 
from  the  series  those  terms  containing  the  lowest  root  of  JJ(nb)  =  0, 
denoted  by  n^.  These  terms  contain  the  transverse-electric  TE^  reso¬ 
nance,  which  is  the  lowest  electromagnetic  waveguide  resonance.  All  other 
self-field  contributions  are  neglected.  This  procedure  is  equivalent  to 
treating  the  beam  as  a  weak  perturbation  of  the  vacuum  TE^  eigenmode. 


We  now  proceed  to  solve  Eq.  (3)  by  writing  £  in  the  form 


.  r  _  +  Ue-iwt  ,  in+iwt 
?  =  e  +  5.  e 


(7) 


where  +,  -  refer  to  "forward"  and  "backward"  waves  respectively.  Substi¬ 
tuting  this  into  Eq.  (3),  we  obtain  the  coupled  equations 


2  12  afi  + 

-n+  +  -j  nz  -  n+si0  +  2 - "TT?. 

w  -  nji  -  l  /R  J 


*°  •  <s-> 


-(n.  -  "«z)2  *  7  +  («+  -  "«z)»e 


a(n+  -  nfi  J 
+  ~1 - T 


-  n^  -  (i  +  m)  /R  J 


-m-t  -’z'-'t 


=  0  , 


(8b) 


+  7  nz  +  +  "2 - 2 


aQ 


— ? — 7 
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0  ,  (8c) 
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l' 


I* 


p 

where  =  w  +  iSi^,  a  =  vl^/y  for  n  =  n^  (o  -  8.4v/yb  ).  Note  that  this 

is  a  closed  system  of  equations,  in  which  only  mode-numbers  t  and  -m  -  l 

appear.  This  is  a  consequence  of  the  complex-conjugate  sign  appearing  in 

Eq.  (3).  Further,  note  that  5*  is  coupled  only  to  c”  ,  and  cj  is  coupled 

only  to  £+  . .  Thus  the  4x4  matrix  of  Eq.  (8)  splits  into  two  2x2 

-m-Jt 

matrices.  The  determinant  of  one  2x2  matrix  gives  the  dispersion  relation 


n 


2  1 
+  “  7 


+  flfle 


u> 


~2 — ? 

1  tvr 


2 


7°z 


(fl+  -  n«z)fi0  -  ~y 


«(n+  -  2V 


7 

’ll 


(l  +  m)2/R2 


2.4  n 
-  v  «2  =  0 


(9) 


The  determinant  of  the  other  2x2  matrix  can  be  obtained  from  Eq.  (9)  by 
letting  w  ♦  -w.  Having  obtained  a  root  for  o>  from  Eq.  (9),  the  correspond¬ 
ing  normal  mode  for  £  can  be  obtained  from  Eqs.  (8a,  8b),  and  is  a  linear 
combination  of  two  terms  with  space-time  dependences  exp(Ue-iwt)  and 
exp[-i(m  +t)e  +  iut]  respectively. 


Multiplying  out  the  dispersion  relation  in  Eq.  (9),  we  obtain  an  eighth 

degree  polynomial  in  w.  Numerical  solution  for  the  sample  parameters  in 

Table  I  gives  the  set  of  curves  shown  in  Fig.  2.  The  parameter  e  in 

Table  I  is  related  to  p  by  u  =  enfl0/2n  .  The  curves  can  be  classified  as 

0  z 

follows: 

2  2  2  2  2  2  2 

Electromagnetic  modes:  w  *  n^  +  l  /R  ,  n^  +  (m  +  t)  /R  , 

Cyclotron  modes:  w  *  (m  +  i)n  +  nQ  ,  to  -  na  , 

Low  frequency  transverse  modes:  w  *>  (m  +  -  u)D  ,  £0  +  u>D  ,  (10) 

ZB  ZB 

2 

where  -  1/4  nz/fi0.  These  approximate  forms  are  based  on  the  assumption 

Q.  »fl„,  which  is  satisfied  for  typical  stellatron  parameters. 

“  z 


t 
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TABLE  I.  HIGH-CURRENT  STELLATRON  PARAMETERS  USED  IN  PARTICLE  SIMULATIONS 


1  m 

8.4  -  9.5  cm 

1  m 

2  cm 

300  A  -  10  kA 
7 

5  kG 
118  G 
0.7 
14-30 

r 


r 


Torus  Major  Radius 
Torus  Minor  Radius 
Beam  Major  Radius 
Beam  Minor  Radius 
Beam  Current  (I) 

Beam  Energy  (y) 

Toroidal  Magnetic  Field 
Vertical  Magnetic  Field 
Quadrupole  Field  (e) 
Quadrupole  Mode-Number 
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r{cm 
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It  can  be  shown  from  Eq.  (9),  and  Is  apparent  from  Figure  2,  that  the 
dispersion  relation  is  symmetrical  about  the  line  i  -  -m/2  (m  is  an  even 
integer,  equal  to  twice  the  number  of  minor  turns  the  helical  conductors 
make  in  going  one  major  turn  around  the  torus).  More  precisely,  if  we 
define  i'  -  £  +  m/2,  then  to  -*•  -to*  as  £*  -*•  -£'.  In  the  following,  we  will 
look  only  at  £'  >  0.  Results  for  £'  <  0  follow  by  symmetry. 

Instability  occurs  due  to  the  intersection  of  the  curve  to  -  (m  +  £)flz 
tog  with  the  electromagnetic  mode  to  •  +  The  instability  can  be 

thought  of  as  a  parametric  process13  in  which  the  quadrupole  field  plays 
the  role  of  a  pump  wave.  If  we  denote  the  three  interacting  modes  by  sub¬ 
scripts  1,2,3,  then  we  get  the  following  frequencies  and  mode-numbers: 

Pump  wave:  toj  =  0  ,  k^  =  m/R  , 

~~2  — 5  ? 

Electromagnetic  wave:  to^  *  +  £  / R  ,  V.^  ~  £/R  , 

Slow  transverse  wave:  •  (m  +  £)fiz  ,  kj  *  (m  +  £ )/R  ,  (11 

In  terms  of  this  picture,  instability  occurs  when  the  usual  matching 
conditions  are  approximately  satisfied  (u>j  +  102  =  wg,  k^  +  k^  *  k^). 

Energy  for  the  instability  is  provided  by  the  beam  via  the  negative-energy 

slow  transverse  mode. 


III.  DERIVATION  OF  INSTABILITY  CRITERION 


An  approximate  expression  for  the  growth  rate  can  be  obtained  by 

2  22  1/2 

letting  w  =  co..  +  6,  where  w..  =  (nf,  +  l  /R  )  '  in  Eq.  (9),  and  keeping 


terms  to  order  6  .  This  gives  the  quadratic  expression 


-2u.11(u,n  -  «.n2)C^11  -  +  fie)V‘ 


2a>2^(u)jj  -  -  £fiz  +  ftg)[  (t  +  m)Oz  - 


+  (u>11  -  £fiz)2  a -  2u11w2fiz 


6  -  o(a>u  -  tflz)' 


x  [ (m  +  i)ny  -  a)n]ne  =  0  . 


(12) 


Defining  the  quantities  y^»  Y2»  Y3  by 

t 1  =  ^liC^n  -  *fi2)(»n  -  +  fig)[(£  +  m)nz  -  0)n]ne  , 

y2  =  a(a)n  -  *fiz)2n0  ,  (13) 

o  2  4 
y3  =  2u,n  U  fiz  , 

the  condition  for  instability  can  be  written 

(Yl  +  y2  -  y3)2  <  4yjY2  ,  (14a) 

or,  equivalently, 

(^1  -  y2  -  y3)2  <  4y2y3  ,  (14b) 

Equations  (12) -( 14)  allow  us  to  make  some  more  exact  statements  about  the 
conditions  for  instability.  From  Eqs.  (13)  and  (14b),  we  see  that  the 
instability  disappears  for  u  =  0,  as  one  would  expect.  Equation  (14a) 
shows  that  (£  +  m)nz  >  (i.e.  Yj  >  0)  is  required  for  instability 

(assuming  >  0;  the  case  <  0  is  discussed  later).  This  means  that 
the  instability  turns  on  when  the  frequency  of  the  transverse  beam  mode  is 
somewhat  above  the  TE^  eigenfrequency ,  rather  than  exactly  equal  to  it. 


To  see  how  the  instability  turns  off  as  £  increases,  we  rewrite  Eq.  (14a) 
to  obtain 


as  the  criterion  for  instability.  The  term  Yj  is  sensitive  to  the  value  of 
£  through  the  factor  (£  +  m)nz  -  As  £  increases,  Yj  eventually 

becomes  large  enough  to  violate  this  inequality,  and  the  instability  dis¬ 
appears  (cf.  Figure  2).  Equation  (15)  also  shows  that  as  v  +  0,  (so  that 
Y2  *  0) ,  and  u  *  0,  the  unstable  region  in  t-space  becomes  narrower.  If  £ 
were  a  continuous  variable,  then  the  instability  would  persist  near 
(£  +  m)£2z  =  Ujj  as  long  as  v  and  u  were  finite.  Since  £  is  discrete  in  a 
toroidal  system,  however,  the  unstable  region  can  fall  between  two  integer 
values  of  £,  and  no  instability  would  be  seen. 


Finally,  from  Eqs.  (12)  and  (14a),  the  peak  growth  rate  as  a  function 
of  u  is  found  when  u  is  chosen  so  that  Ti  +  Y2-Y3s0.  The  growth  rate 
r  is  then  given  by 


r2  -  a 
r  ‘  1 


[  (m  +  £)n  -  -  «J_) 


“11^11 


For  the  parameters  of  Fig.  2,  this  equation  predicts  r  =  1.06  x  10"2  for 
£  =  0,  in  good  agreement  with  the  exact  result  of  1.12  x  10"?  which 
occurs  for  e  =  0.8  (y  -  380). 


Thus  far,  we  have  assumed  >  0.  If  £la  is  negative,  then  Eq.  (14b) 
shows  that  the  instability  disappears.  This  is  to  be  expected  since  the 
slow  transverse  wave  then  becomes  a  fast,  positive  energy  wave.  At  the 
same  time,  the  cyclotron  wave  oj  ■  (m  +  £)nz  +  becomes  a  slow  wave,  with 
the  potential  for  an  unstable  interaction.  We  have  not  examined  this  case, 
however,  since  previous  calculations3  have  shown  that  the  parameter  space 
in  which  single  particle  orbits  are  stable  shrinks  greatly  when  the  sign  of 
is  opposite  to  that  of  m. 


IV.  COMPARISON  WITH  NUMERICAL  SIMULATIONS 


To  verify  the  above  analytic  calculations,  we  have  performed  three- 
dimensional  numerical  particle  simulations  with  the  code  IVORY.  For  these 
simulations,  field  quantities  are  assumed  to  have  the  form 

F(  r,e  ,z,t)  =  l  fn(r,z,t)  e1n6  .  (17) 

n 

In  the  r-z  plane,  field  quantities  are  represented  on  a  two-dimensional 
spatial  mesh.  The  self-consistent  fields  of  the  beam  are  advanced  in  time 
using  the  full  Maxwell's  equations.  Particles  are  advanced  using  the  full 
Lorentz  force  equations.  The  stellatron  fields  are  computed  from  analytic 
expressions.3  In  deciding  on  parameters  for  the  simulations,  computing 
costs  constrain  us  to  choose  cases  which  minimize  the  running  time  and 
storage  requirements.  Thus,  we  concentrate  on  cases  with  large  expected 
growth  rates.  In  addition,  we  choose  either  4  =  0,  or  choose  m  +  4  to  be 
a  small  integer  multiple  of  4.  This  minimizes  the  number  of  particles 
needed  to  resolve  the  different  mode-numbers.  For  4*0,  for  example,  we 
see  from  Eqs.  (11)  and  (17)  that  only  mode-numbers  0,  +m,  must  be  repre¬ 
sented  in  the  simulations. 

Choosing  4  =  0,  we  used  the  parameters  in  Table  I  with  b  =  8.8  cm  and  a 
10  kA  beam  current,  and  performed  simulations  for  different  values  of  the 
quadrupole  mode-number  m.  The  growth-rates  obtained  are  plotted  in  Fig.  3 
versus  those  obtained  from  Eq.  (9).  We  also  show  the  growth  rate  from  a 
simulation  for  m  =  22  in  which  the  beam  was  represented  by  a  string  of 
rigid  disks  in  the  code,  instead  of  the  more  realistic  particle  representa¬ 
tion.  There  are  at  least  two  possible  reasons  why  agreement  with  the 
analytic  result  is  better  for  the  rigid  disk  simulation.  First,  the  ini¬ 
tial  field  energy  level  is  lower  for  the  latter  due  to  the  absence  of 
internal  degrees  of  freedom  of  the  beam,  so  that  there  is  a  longer  visible 
period  of  linear  growth.  In  the  simulations  with  particles,  the  field 
energy  only  increases  only  about  one  order  of  magnitude  over  its  initial 
value  before  saturating,  as  shown  in  Figure  4.  This  means  that  the  growth 
measured  may  not  be  truly  exponential.  Second,  in  the  particle  simulations 
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Figure  4.  Energy  inn  =  22  field  components  (n  =  toroidal  mode- 
number)  versus  time  showing  development  of  l  =  0, 
m  =22  instability  in  Figure  2.  This  plot  was  generated 
by  the  simulation-code  IVORY. 


there  is  a  spread  In  the  energy  of  the  particles  due  to  space-charge, 
giving  rise  to  a  spread  in  transverse  oscillation  frequencies.  This  may 
have  a  stabilizing  effect. 

There  is  good  agreement  between  simulations  and  theory  as  to  the 
turn-on  and  turn-off  of  the  instability.  The  TE^  cutoff  frequency  for  the 
drift-tube  is  *  1.84/b  -  0.21.  Figure  3  shows  that,  in  agreement 
with  the  analytic  prediction  (Sec.  Ill)  Instability  sets  in  for 
m  >  n^/G  *  21.  No  instability  is  seen  during  the  length  of  the  simula¬ 
tions  for  m  28  (we  ran  cases  for  m  *  28,  30). 

Further  evidence  of  the  electromagnetic  character  of  the  instability 
comes  from  two  simulations  where  the  minor  cross-section  of  the  drift -tube 
was  varied  keeping  m  fixed.  For  m  «  20,  I  =  10  kA,  b  =  8.8  cm  (other 
parameters  as  in  Table  I),  no  growth  was  observed,  since  m/R  <  1.84/b.  On 
increasing  the  minor  radius  to  9.5  cm,  however,  strong  growth  was  observed. 
In  addition,  a  field  contour  plot  of  the  n  *  0  component  of  the  perturbed 
B0,  shown  in  Figure  5,  reveals  a  TE^  structure. 

In  order  to  see  whether  the  aspect  ratio  of  the  torus  affects  the 
instability,  we  increased  the  major  radius  R  first  to  10  m  and  then  to 
100  m,  keeping  m/R  fixed  at  0.22  and  choosing  i  =  0.  No  significant  change 
in  the  growth-rate  was  observed  in  the  simulations.  This  is  in  contrast  to 
the  negativp-mass  instability,14  where  the  growth  rate  typically  falls  off 
as  1/R. 

1/2 

Next,  we  tested  the  v  dependence  of  the  growth  rate  on  the  beam 
current  predicted  by  Eq.  (16).  The  parameters  in  Table  I  with  b  =  9.0  cm, 
m  =  22  were  used  to  perform  simulations  at  300  A,  1  kA  and  10  kA  beam 
currents. The  results  in  Figure  6  show  that  the  predicted  scaling  is 
supported  by  the  simulations. 


1.297  x  10 


Thus  far,  we  have  reported  on  results  for  i  =  0.  Now,  we  look  at  a 
case  for  which  I  -  m  *  14.  The  dispersion  diagram  for  this  case,  obtained 
from  Eq.  (9),  is  similar  to  that  in  Figure  2,  except  that  one  of  the 
unstable  regions  occurs  around  I  =  14.  The  simulation  parameters  used  are 
those  in  Table  I,  with  I  =  10  kA,  b  =  8.4  cm  chosen  to  enhance  the  growth 
rate.  From  Eqs.  (11),  (17),  we  see  that  mode-numbers  0,  ±z,  ±(i  +  m)  must 
be  represented  in  the  simulation  code  {i  =  0  must  always  be  present  to 
represent  the  equilibrium  fields).  Field  energy  in  the  modes  n|  =  t, 

(i  +  m)  is  observed  to  grow  exponentially  at  a  rate  r  =  6.3  x  10"3, 
compared  to  the  analytic  rate  7.8  x  10“^.  In  this  case  it  is  the  |n|  = 

14  field  plots  which  show  a  TE  mode  character  similar  to  that  in  Figure  4. 


The  nonlinear  development  of  the  instability  is  an  important  issue.  To 
address  it  completely,  many  modes  would  have  to  be  kept  in  the  simulation 
code,  since  nonlinear  effects  give  rise  to  the  generation  of  modes  other 
than  those  involved  in  the  linear  growth  stage.  However,  we  believe  that 
the  simulation  results  with  just  the  linear  modes  present  may  give  a  good 
guide  to  the  nonlinear  development  of  the  instability  for  the  following 
reasons.  First,  the  quadrupole  field  gives  a  initial  perturbation  to  the 
mode  number  m,  so  that  the  fields  in  this  mode  are  not  growing  from  random 
noise  (see  Figure  4).  Second,  although  wave-wave  interactions  are  not 
treated  correctly  if  we  do  not  include  other  modes,  the  wave-particle 
interactions  are  treated  nonlinearly.  The  simulations  show  that  in  those 
cases  with  large  linear  growth  rates,  which  we  can  therefore  afford  to  run 
to  "saturation",  the  wave  amplitude  grows  until  the  beam  strikes  the  wall 
(see  Figure  4).  As  a  result,  most  of  the  beam  particles  are  lost,  leaving 
a  large  amplitude  TE^  wave  in  the  drift-tube. 

Finally,  we  look  at  what  our  analysis  predicts  for  the  stellatron 
experiment  presently  under  way  at  the  University  of  California  at  Irvine 
(UCI).15*16  The  parameters  for  this  experiment  are  tabulated  in  Table 
II.  Because  of  the  low  current,  the  width  of  the  instability  is  quite 
narrow,  extending  over  at  most  one  z-number.  As  a  result,  small  changes  in 
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SIMULATIONS 


Figure  6.  Growth  rate  of  instability  versus  beam  current  for 
Table  I  parameters  with  m  =  22,  b  =  9.0.  Simulation 
results  are  compared  to  theoretical  prediction. 


V.  SUMMARY  AND  CONCLUSIONS 


We  have  derived  a  dispersion  relation  for  a  parametric  electromagnetic 
instability  in  a  stellatron  accelerator.  The  instability  arises  from  the 
interaction  between  the  quadrupole  winding,  a  negative-energy  transverse 
wave  on  the  beam,  and  an  electromagnetic  wave-guide  mode.  The  growth  rate 
of  the  instability  is  independent  of  the  radius  of  the  toroidal  drift- 
tube.  It  therefore  occurs  in  straight  as  well  as  toroidal  systems.  Three- 
dimensional  numerical  simulations  of  the  stellatron  have  been  carried  out 
with  the  code  IVORY.  The  simulated  linear  growth  rates  and  conditions  for 
the  onset  of  the  instability  are  in  reasonably  good  agreement  with  the 
analytic  model.  The  simulations  show  strong  disruption  of  the  beam  in  the 
nonlinear  regime,  leading  to  loss  of  current.  Our  calculations  predict 
some  growth  of  the  instability  for  the  parameters  of  the  UCI  stellatron 
experiment.  The  instability  may  be  a  factor  in  limiting  the  beam  current 
in  this  experiment. 
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